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In the paper, we discuss the studies of mathematical models of 
diffusion scattering of waves in the phase space, and relation of these 
models with quantum mechanics. In the previous works it is shown 
that in these models of classical scattering process of waves, the quan- 
tum mechanical description arises as the asymptotics after a small 
time. In this respect, the proposed models can be considered as ex- 
amples in which the quantum descriptions arise as approximate ones 
for certain hypothetical reality. The deviation between the proposed 
models and the quantum ones can arise, for example, for processes 
with rapidly changing potential function. Under its action the diffu- 
sion scattering process of waves will go out from the states described 
by quantum mechanics. 

In the paper it is shown that the proposed models of diffusion scat- 
tering of waves possess the property of gauge invariance. This implies 
that they are described similarly in all inertial coordinate systems, 
i. e., they are invariant under the Galileo transformations. 

We propose a program of further research. 

1 Introduction 



Usually description of quantum systems is constructed by using formal quan- 
tization procedures, based on the classical description of the corresponding 
mechanical systems. The search for the sense of these procedures attracted 
many physicists, starting with A. Einstein. 

The interest to this subject periodically decreased and revived again. In 
this direction, one can mention the von Neumann theorem, proved during the 
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period of formation of quantum mechanics, on impossibility of description of 
quantum mechanics by introducing hidden parameters [T]. Nevertheless, in 
the 50s, in the works of D. Bohm and L. de Broglie [2], one proposed a 
model of quantum mechanics with hidden parameters, not satisfying some 
conditions of the von Neumann theorem and possessing a strange property 
of far-action. In the 60-s, in the work of E. Nelson [3], one proposed a prob- 
abilistic approach given the name "stochastic quantization". The subject 
of foundation of quantum mechanics excited many specialists (for example, 
D. I. Blohintsev [1], V. P. Maslov [5], K. Popper etc.), and it was dis- 
cussed in their publications. There are works in which one makes a detailed 
analysis of the problem of introducing hidden parameters into quantum me- 
chanics. They include a widely known work of Bell [7j on introducing hidden 
parameters and non-locality of quantum mechanics. An interesting analysis 
of this work is given in [8]. The paper |9] contains a huge list of literature 
on foundations of quantum mechanics, and provides a classification of these 
works. The history of the discussion around the subject "foundation of quan- 
tum mechanics" and attitude to this subject of "traditional physicists" are 
remarkably described in the book by K. Popper [6]. 

The small popularity of alternative approaches to foundation of quantum 
mechanics among the working physicists is usually related to the fact that 
they did not yet give serious new results. They also did not give more con- 
venience in computations and heuristics. However, recently the alternative 
approaches cause again an intent attention related to the problems and possi- 
bilities of quantum optics, as well as the problem of construction of quantum 
computers. 

In the present paper we discuss the research on construction of models of 
diffusion scattering of waves in phase space. I have been studying this subject 
during last years (TUl [HI [121 [13] ■ In these models the quantum description 
of processes arises as an approximate one, asymptotical for large values of 
certain coefficients of the model. 

In the papers mentioned above one makes an attempt to construct a model 
of quantum observables on the base of wave functions on the phase space. 
Note that in quantum mechanics, the wave function depends either only on 
coordinates or only on momenta, while in the present approach one considers 
wave functions depending both on coordinates and on momenta. This model 
is based also on the following observation. In quantum mechanics, the phase 
of the wave function of a particle (the natural hidden parameter) changes 
in time even for stationary states with very high velocity (if one takes into 
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account the stationary energy). This velocity is such that a transfer of the 
particle with even small (non-relativistic) velocities can cause considerable 
changes in the phase of wave function because of the relativistic effect of more 
slow inner processes of a moving particle. Already taking into account this 
effect leads to non-commutativity of the action of coordinate and momentum 
shifts on the wave function. Note once again that in the proposed model 
one considers wave functions on the phase (and not configuration) space, 
and one assumes that the particle is in a diffusion process causing random 
shifts of the wave both by coordinates and by momenta. It is shown that the 
classical model of scattering of the wave, taking into account the assumptions 
described above, yields to arising quantum effects. 

In the further sections of the present paper we speak in more detail on 
the obtained results and directions of further research. 

2 The results obtained earlier 

In the paper [10] one introduces some assumptions on the process of obser- 
vation of quantum phenomena, including introduction of hidden parameters, 
action of the group of motion in the region of hidden parameters, and aver- 
aging observations due to small random (diffusion) motions of the observed 
object. By an observable we mean, as in classical mechanics, an arbitrary 
integrable function f{x,p) on the phase space {x,p) G -R^", where x is the 
coordinate, p is the momentum. If p{x,p) is the density of probability distri- 
bution of the position of a particle in the phase space, then the mathematical 
expectation (mean) / of an observable / is given by the standard formula: 



Below it is assumed that in experiments, not all distributions p{x,p) are 
realized, but only those of the form p{x,p) = |<^(x,p)p, where (p{x,p) is a 
wave function averaged in the diffusion process, given in the form of a complex 
valued function ip{x,p) on the phase space. It is shown that the functions 
of the form (p{x,p), form a linear subspace T-L of "stationary" (averaged) 
wave functions in the space of all square integrable functions on i?^". Since 
= (p*(p^ where (p* is the complex conjugate function to the function 
ip, the mean value of the observable / on averaged densities of probability 
distribution yields the following quadratic form on the space of G 
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f= j f{x,p)p{x,p)dxdp= J f{x,p)(p*{x,p) <^{x,p)dxdp = {^,Af(p), 

^2n ^2n 



where by Aj we denote the hnear operator on the space "H giving this 
quadratic form. 

The introduced operator Af is called the operator of the observable /. It 
is natural that the spectrum of this linear operator corresponds to possible 
values of observations for the observable / under the assumptions made. 

In the paper [10] we have found the expression for this operator for any 
observable (function of coordinates and momenta), depending on the ratio 
a/6 of diffusion coefficients, with respect to coordinates a and momenta b, 
of the averaging process of wave functions. It is shown that the usual linear 
operator of a quantum observable does not coincide with the one constructed 
in the paper, but differs by smoothing of the potential energy function with 
respect to the normal distribution with the normal deviation equal to ha /2b, 
where h is the Planck constant. Assuming that this difference yields the shift 
of the spectrum of the hydrogen atom observed in the Lamb experiment, we 
give an estimate of the quantity a/b. 

A big advantage of the considered approach is also the possibility to 
express, for each wave function of the system, the corresponding density 
of probability distribution in the phase space p{x,p) = \ip{x,p)\'^. For the 
first time this problem was solved by Wigner [H], but he has constructed 
"quasidistributions" on the phase space which can be negative and hence have 
no physical sense. And here we have a probability density distribution which 
is the result of smoothing of Wigner's "quasidistribution" with respect to 
the normal distribution with the normal deviation equal to ha /2b. Smoothed 
Wigner's distributions were first considered by Husimi [I5], but the sense of 
the smoothing parameters was unclear. 

At the end of the paper yLO] we have posed the problems of generalization 
of the results to the relativistically invariant case, taking into account the spin 
of the particles, more general phase manifolds, and description of dynamics 
of observable quantities. 

The papers |Tll |12l US] are devoted to solution of the latter problem. 

In the papers [HI [12] (in [11] the results of [12] have been announced), 
continuing the work [10], we consider the classical model of the diffusion 
process for a wave (complex valued) function ip{x,p,t) on the phase space 
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{x,p) G -R^" at the moment of time t. It is assumed that the wave function 
ip{x,p,t) at the moment t satisfies the following differential equation: 



where A,,^ = « E - + & E + — V^, (2) 

where H{x,p) is the Hamilton function; and 6^ are the diffusion coefficients 
with respect to coordinates and momenta respectively, and h is the Planck 
constant. 

The analysis of this equation has shown (see [12], Theorems 4 and 5) that 
in this model the motion splits into rapid and slow ones. After the rapid mo- 
tion, at the time of order h/ {ah), starting from an arbitrary wave function on 
the phase space, the system goes to a function belonging to certain special 
subspace of "stationary" wave functions for the diffusion process. The ele- 
ments of this subspace are parameterized by wave functions depending only 
on coordinates. The slow motion takes place already in this subspace and is 
described by the Schrodinger equation, in which in the right hand side we 
have the operator coinciding with the usual quantum mechanics Hamilton 
operator up to summands of order ah/h. 

Thus, already in these papers it is shown that the quantum mechani- 
cal description of processes can arise as the approximate description of the 
classical diffusion of waves in the phase space. For the model considered in 
the paper, this approximation arises when the Hamilton function has a small 
change with a change of coordinates, momenta and time in intervals of length 
of order defined by the Planck constant and diffusion intensities. 

Assuming the heat reason of the diffusions, in the paper we estimate 
the diffusion coefficients and the transition time h/{ab) from the classical 
description of the process in which the Heisenberg indeterminacy principle 
in general does not hold, to the quantum description in which the Heisenberg 
principle already holds. The transition time has order 1/T ■ 10~^^ sec, where 
T is the temperature of the medium. 

Another interesting result of the paper [11] is that the solution of equation 
(II]) can be represented as a path integral, but not with respect to the Feynman 
"measure" [16], whose sense is mathematically not very much clear, but with 
respect to the probability measure (analogous to the Wiener measure) for 
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the Brownian motion given by the Fokker-Planck equation of the form 



Here f{x,p,t) is the probabihty density of the position of the Brownian 
particle in the phase space at the moment of time t. In this case, the sense 
of the path integral can be better substantiated. 

Generalization of equation ([1]) to the relativistic case meets some difficul- 
ties, because of the presence of diffusion with respect to coordinates in this 
model. Such diffusions assume unbounded velocity in the diffusion jumps. 
Hence the next step in our investigations was construction of a model of 
scattering of waves in the phase space, in which the diffusion takes place 
only with respect to momenta, because of the collision with particles of the 
medium in the heat equilibrium. 

In the paper [13], instead of equation (|3]) we consider the Kramers equa- 
tion [T7], [IB] of the form 

^ = y(^^-PiK)+^y^(p.f + kTm^) (4) 

dt ^i\dxj dpj mdxjj J^idpjX ^ dpjj' 

where f{x,p,t) is the probability density of a particle in the phase space at 
the moment of time t; m is the mass of the particle; V{x) is the potential 
function of external forces acting on the particle; 7 = /3/m is the resistance 
coefficient of the medium in which the particle moves, per unit of its mass; 
k is the Boltzmann constant; T is the temperature of the medium. 

Then, instead of equation ([T]) for the wave function (p{x,p, t), we consider 
the modified Kramers equation of the form 

^ = A^ + 7Sy^, (5) 

B*' = 1 1- ( (ft + a^) + . (7) 

Equation ([5]) is obtained from the Kramers equation (j4]) by adding to the 
right hand side of the summand of the form —i/h{mc^ + V —p^ / {2m))ip, and 
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the replacement, in the diffusion operator, of multiphcation of the function 
V9 by Pj by the action of the operator {pj + ihd/dxj) on the function ip. 

Adding the summand —i/h{mc^ + V — p^ / {2m))ip is related with the 
additional physical requirement that the wave function at the point {x,p) 
oscillates harmonically with frequency l/h{mc^ + V — p^ /{2m)) in time. 

The requirement of harmonic oscillating of the wave function at the 
point with the large frequency l/h{mc^ + V — p^/(2m)), in the case 

when mc^ is much greater than V, leads to the fact that the shift of the wave 
function with respect to the coordinate Xj with conservation of the proper 
time at the point (x, p) yields the phase shift in the oscillation of the function 
Lf. And the operator of infinitely small shift d/dxj is changed by the operator 
d/dxj —ipj/h. (For a more detailed explanation, see [12] •) Respectively, if 
we multiply this operator by ih, then we obtain the operator pj + ihd/dxj, 
used in the modified diffusion operator B. 

For equation ([5]), in [13] we obtain results similar to that of the paper 
[12]. It is shown that also in this case, the process described by equation (|5]), 
for large 7 = /3/m passes several stages. During the first rapid stage, the 
wave function goes to a "stationary" state of the same form as for equation 
(|5]). At the second, slow stage, the wave function evolves in the subspace 
of "stationary" states subject to the Schrodinger equation. Besides that, 
it is shown that at the third stage, the dissipation of the process leads to 
decoherence of the wave function, and any superposition of states comes to 
one of eigenstates of the Hamilton operator. 

In the paper [13], it is shown also that if, on the contrary, the medium 
resistance per unit of mass of the particle 7 = /3/m is small, and in equation 
(|5]) one can neglect the summand with the factor 7, then in the considered 
model, the density of the probability distribution p = lyjp satisfies the stan- 
dard Liouville equation 



as in classical statistical mechanics. 

3 Gauge transformations 

In this section we introduce and discuss the notion of gauge invariance for 
equation 



dp 
di 
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According to the approach exposed in [13], the density of probabihty 
distribution p{x,p,t) of a quantum particle whose state at the moment of 
time t is given by the wave function (p{x,p,t), is proportional to = 
(p{x,p,t)ip*{x,p,t). This implies that the replacement of a wave function ip 
by the wave function of the form exp{ig/h)ip, where g = g{x,p,t) is an 
arbitrary real valued function, does not change the density of the probability 
distribution p{x,p,t). Such a transform of wave function is usually called a 
gauge transform. 

Let us look how equation ([5]) changes under this gauge transform. To 
this end, let us write out equation ([5]) in a more general form. Let us write 
in it, instead of the differentiation operators d/dpj of the function (p, the 
operator = d/dpj + iBj/h, instead of the operators d/dxj — ipj/h, the 
operator Dj = d/dxj+iAj/h, and instead of the operator d/dt+iH/h, where 
H = mc^ +p'^ /(2m) + V , let us write the operator Dq = d/dt + iAo/h, where 
Aj,Ao,Bj are functions of x,p, and t for j = l,...,n. In these notations, 
equation ([5]) will take the form 

Do^ = E - g^D^^) + 7 E i^nn]^ + kTmD^^) . (9) 

By a gauge transform of equation we call the following transform of 
the function ip and the potentials Aj, Ao, Bj, for j = 1, n: 

(10) 

l,...,n; 

l,...,n. (11) 

It is not difficult to see that after the substitution ( ITOl) into equation 
dHI), replacement (fTTl) . and dividing both parts of the obtained equality by 
exp{—{i/h)g), the form of equation Q) will not change. 

Geometrically, gauge transformation corresponds to transfer to another 
trivialization of a complex line bundle over the phase space, in which a form 
of linear connection is chosen, defining parallel transport of the vectors of 
the bundle along trajectories in the phase space. 



I — y (p' = exp{--g)<p; 



Ao ^ A'q = Ao + 



dg_ 

dg 



A^ ^ 4 = A,- + where j 
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In the particular case for equation (jS]), the potentials read 



Aq = H{x,p) = E + V; Aj = -pf, Bj = for j = 1, n. 

Understanding the physical sense of the potentials in the general case 
for equation ([9]), requires separate investigation. For the Dirac equation, 
potentials of gauge invariance are usually related with the potentials of elec- 
tromagnetic field. 

4 The Galileo invariance 

In this section we study the change of equation ([5]) under the transfer to a 
coordinate system moving uniformly with respect to the initial coordinate 
system, with the velocity u. The diffusion equation (jl]) is not invariant with 
respect to Galileo transforms under transfer to new inertial coordinate system 
moving with constant velocity u with respect to the old one. 

The aim of this section is to study invariance of equation ([5]) for a free 
particle {V = 0) with respect to Galileo transforms, with gauge transforms 
of the wave function. 

By definition of Galileo transforms, the new coordinate system is ex- 
pressed through the old one by the following formulas: 

t' = t; x' = X — ut; p' = p — mu; 

r-,/ P'"^ ip — muY p^ rnu^ ^ mu^ 

E = — = = pu + — - = E~pu+ — -. 12 

2m 2m 2m 2 2 

Respectively, the old coordinates are expressed through the new ones by 
the following formulas: 

t = t'] X = x' + ut] p = p' + mu] 

E=^=^-^±^ = ^ + p'u + ^ = E'+p'u + ^. (13) 
2m 2m 2m 2 ^ 2 ^ ' 

Substituting these expressions into equation (jSj), with the use of relations 
([6]) and ([7]), we obtain: 

dip ^ dip _ ^/dV dip p'j + muj/ d ,p'j+muj\ \ 
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-j^\E'+p'u+ — + V^^ 

whence, after simple algebraic transformations, we obtain: 

(9v? _ ^/dVdip p'j / d ,p'j + muj\ 
~{\dx'j dp'j m\dx'j ^ h / 

If in the obtained equation one makes the substitution ip = exp{{i/h)g)ip' , 
where g = mux' + muH' /2, then (after the gauge transform) we obtain the 
equation 



d^' _ ^/dV dif' _p'j / d _.p'j 
9t' ~{\ dx'j dp'j m \ dx'j Ti 



-\{E' + V)v' 

which coincides with equation (E]). Thus, we have proved the Galileo invari- 
ance of equation ([5]). 



5 Program of further research 

In this section we list directions of further research and sketch approaches to 
the stated problems. 

5.1. Comparison of the model of scattering of waves 
with the quantum model 

In order to compare exactness of the model described by equation (|5]), 
with the standard quantum mechanical model, one should find the situation 
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in which these models give essentially different results. Such a situation can 
arise, for example, if one considers the process with the rapidly changing in 
time potential function V{x, t). Such a potential can prevent a wave function 
of equation ([5]) from transfer, during the time of the transition process, to 
the "stationary" one. As a result, a solution of equation ([5]) can differ from 
a solution of the Schrodinger equation. 

In order to check this, consider, for example, the potential function V = 
Vo(x) + Vi{x) cos(wt) for u — > oo. 

Mechanical and quantum mechanical systems with such potential were 
studied in many papers, for example, [HI |20l |2ll |22l [231 [211 [25] . The physical 
problem in which such a quantum model arises, is a charged particle in 
external force field and in a laser row. 

Equation ([5|) with this potential reads 

d(p 



where 



Aif + ^Bif, (14) 



fd{Vo + Vi cos{ut)) dip Pj dip 
\ dxj dpj mdxj 



j=i ^ ^-^j ^yj ^-^j 

n2 



j^y'"^ I .-u I ^icos(wt) - (15) 



mc^ + Vn 



and B^ = |:|-((p,+«,|-)^ + .r™g. (16) 

One should study solutions of this equation for large u and compare these 
solutions with solutions of the quantum system. 

5.2. The study of the scattering process of mixed 
waves and computation of the time of the transition 
process to stationary mixed state of heat equiUbrium 

Another problem which one would like to study is the behavior of the 
process for mixed waves of the form ip{x,p, t, ^), where ^ e D is an additional 
parameter, and the distribution p{x, p, t) in the phase space at the moment 
t for a particle whose state is described by a wave function ip{x,p,t,C,), is 
proportional to the function \ip{x,p,t,^)\'^d^, i. e. 

Id P, t, ^)ip*{x,p, t, ^)d^ 



p{x,p,t) 



Sa^r. Jd <^{x, p, t, p, t, Cjd^dxdp' 
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Also here one assumes that the evolution of the wave function in time goes 
according to equation for each fixed ^ E D. 

Another equivalent way to describe this process, familiar in quantum me- 
chanics, is to consider self-adjoint operators p on functions on the phase 
space i?^" with the kernel of the operator of the form p{x,p] x' ,p' ,t) = 
J^ip{x,p,t,C,)ip*{x',p',t,^)dC,. Note that any positive self-adjoint operator 
p on the space of functions can be reduced to diagonal form and there- 
fore to the form above. Positive self-adjoint operators with trace unity are 
called operators of density of states. Then the density of probability distri- 
bution p{x,p,t) = p{x,p; x,p,t)/Trp, where Trp = Jj^2,i p{x,p; x,p,t)dxdp is 
the trace of the operator p. The evolution of the operator of density of state 
p in time is given by the equation 

^ = Vp + pV*- pTr{Vp + pV*), 

where V is the operator expressed by the right hand side of equation ([5]), 
and V* is the adjoint operator. Expression with the trace Tr stands in this 
equation in order to make the trace of operator of density of state p equal to 
one at each moment of time. 

This is a nonlinear equation. One should investigate whether it has a 
unique stationary state, determine the form of this stationary state (the 
state of heat equilibrium), and estimate the time of transition process to this 
stationary state. 

5.3. Generalization of the model with account of spin 
of a particle and the requirement of relativistic invari- 
ance 

Let M = i?^ be the Minkowsky space-time, P = he the space of 
momenta, and B = M x P he the phase space-time, on which the Lorentz 
group naturally acts (if one fixes the stationary mass m of the particle). 
The same space has an action of the commutative group of coordinate shifts 
preserving the proper time at each point of the phase space, and of the one- 
parameter group of shifts of proper time at each point of the phase space. 
Together these groups define an action of the Poincare group P on the space 
B. 

In this new model, we propose to consider the values of the wave function 
(p not in the field of complex numbers C, but in certain Euclidean vector 
space F over the field of complex numbers, with an action by unitary linear 
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operators of the group SU{2, C), the two- fold covering of the rotation group 
5*0(3) of three-dimensional space, acting on the phase space R^. The proba- 
bihty distribution p{x,p,t) of position of a particle in the phase space at the 
moment of time t is again assumed to be proportional to \ip{x,p,t)\'^. 

The group SU{2, C) is a subgroup in the group SL{2, C), where SL{2, C) 
is the group of two-dimensional complex matrices with determinant equal to 
1. The group SL{2, C) is the two-fold covering of the Lorentz group L. Thus, 
we have a commutative diagram of homomorphisms of groups: 



where P is the two-fold covering group for the Poincare group. 

Further one considers the bundle pr : F x B ^ B with fiber F over the 
phase space-time B. The Poincare group P acts on the base B. The action 
of its subgroup 5*0(3) G L on B by rotations with respect to the coordinate 
origin lifts to the compatible action of the group SU{2,C) in the fiber F 
over the origin point in B. Then, this action can be uniquely extended to an 
action of the group P on the bundle F x B, compatible with the action of the 
group P on the base B. The compatibility of the actions of the groups on 
the bundle means that for any g & P, the following diagram is commutative: 



is commutative. 

The uniqueness of the lift of the action of the Poincare group from B to 
the action of the group P on the bundle F x B is understood up to a choice 
of trivialization of this bundle. 

If a, 6 e S are two points of the base (the phase space-time), then one 
uniquely defines an element ha,b € -P of the Poincare group, of the parallel 



SU{2,C) C SL{2,C) C P 
S0{3) C L CP, 




B H B. 



Here also, ii g E SU (2, C) C P, then the diagram 



F X C F X B 
F X c F xB. 
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transport of the coordinate system from the point a to the point b. The 
action of the element ha^b hfts uniquely to the action of an element ha^b ^ P 
on the bundle F x B. This action transfers elements of the fiber F over a to 
elements of the fiber over b. Let us call this action by the parallel transport 
of elements of the fiber along the vector ab. Further, this definition allows us 
to define the parallel transport in the bundle F x B along any curve in the 
base. 

In the considered model, the wave function ip at the moment of time t 
is given by a function on the phase space of the form ip : ^ F. The 
evolution of the wave function in time is defined by the condition that it is 
simultaneously in several motions: 

1) The vector ip{x,p) G F is parallel transported along the trajectory 
in the phase space; the trajectory is defined by a random Brownian process 
according to certain diffusion equation, for example, the Kramers equation. 

2) The vector (p at each point {x,p), in the coordinate system related to 
this point, rotates with the constant angular velocity u = mc^/h in the fiber 
F over this point in the proper time related to this point; the direction of 
the rotation axis J^^p G su{2,C), in the stationary (laboratory) coordinate 
system, transforms from one point to another in the same way as the direction 
of the angular momentum. 

The value of the wave function at the point (x, p) at the moment {t + At) 
is defined by the mean value of the vectors (p over all trajectories ending at 
the point {x,p) of the phase space at the moment (t + At). 

One should construct the differential equation corresponding to this model, 
and study it. 

5.4. Scattering of waves on the phase space and in- 
teraction with electromagnetic field 

This problem is related to introducing interaction with electromagnetic 
field into the model. Such introducing could be made by analogy to its in- 
troducing into the Dirac equation. As it was shown in equation (|9]), on this 
way potentials arise depending also on the momentum, in contrast with the 
vector potential of the electromagnetic field which depends only on coordi- 
nates and time. Determining the sense of vector potentials depending on 
momenta, also requires a separate investigation. 

Acknowledgements: I am grateful to professor A. V. Stoyanovsky, who 
translated this paper to English. 
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